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Abstract 

An analysis of extension of Hamiltonian operators from lower order to higher order 
of matrix paves a way for constructing Hamiltonian pairs which may result in hered- 
itary operators. Based on a specific choice of Hamiltonian operators of lower order, 
new local bi-Hamiltonian coupled KdV systems are proposed. As a consequence of 
bi-Hamiltonian structure, they all possess infinitely many symmetries and infinitely 
many conserved densities. 



1 Introduction 

Bi-Hamiltonian structure has been identified as one of the basic mechanisms supplying con- 
servation laws and commuting flows and even constructing multi-gap solutions for a given 
nonlinear system. From differential geometric point of view, bi-Hamiltonian structure can 
lead to an interesting application which yields a key to describe integrable surfaces embed- 
ded in pseudo-Euclidean spaces. For example, in many cases local bi-Hamiltonian structure 
provides additional information about some nonlocal Hamiltonian structures closely con- 
nected with metrics of constant curvatures. There exist other applications in handling 
quantization of Poisson brackets, and studying the problem of instability, for instance, in 
fluid mechanics. What we would like to develop in this paper is to establish more local 
bi-Hamiltonian structures which can engender hierarchies of coupled KdV systems. 

Let u be a dependent variable u = (u l , ■ ■ • , u q ) T , where u l , 1 < i < q, depend on the 
spatial variable x = (x±, ■ ■ ■ ,x p ) T and on the temporal variable t. We use A r (r > 1) to 
denote the space of r dimensional column function vectors depending on u itself and its 
derivatives with respect to the spatial variable x. Sometimes we write A r {u) in order to 
specify the dependent variable u. 



A natural inner product over A r (u) is given by 



<a,(3>= l^ p (a) T (3dx, a,p € A r . (1.1) 
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Moreover on A q {u) (note that the dimensions of both a function vector in A q (u) and the 
dependent variable u are the same), a Lie product can be defined as follows 

[K, S] = ^(K(u + eS) - S(u + e*0)| e _ 0> K, S e A q . (1.2) 

Two function vectors K,S £ A q (u) are called to be commutative if [K, S] = 0. A funda- 
mental conception that we need is Hamiltonian operators, which is shown in the following 
definition. 

Definition 1.1 A linear skew-symmetric operator J{u) : A q {u) — ► A q {u) is called a 
Hamiltonian operator or to be Hamiltonian, if the Jacobi identity 

< a, J'(u)[J(u)P}j > +cycle(a,(3,j) =0 (1.3) 

holds for all a, f3, 7 G A q (u). A pair of operators J(u), M(u) : A q {u) — > A q (u) is called a 
Hamiltonian pair if J{u) + cM{u) is always Hamiltonian for any constant c. 

Associated with a given Hamiltonian operator J(u), the Poisson bracket is defined to 

be 

{H 1 ,H 2 } J = J( 6 -^ fj 6 J^dx, H\,H 2 £ A, (1.4) 

where A consists of functionals H = f Hdx, H £ A. Two functionals H±,H 2 are called 
to be commutative under the Poisson bracket associated with J, if {Hi,H 2 }j = 0. A 
Hamiltonian operator J : A q — ► A q has a nice property 

{ jfiy^ ] = /AilAk,H u H 2 ,A, (1.5) 
ou ou ou 

which gives rise to an important relation between symmetries and conserved densities for 
a Hamiltonian system with the Hamiltonian operator J. 

If we have a Hamiltonian pair J and M, one of which is invertible, for example J, then 
$ := MJ" 1 is a hereditary operator [E] Q. Furthermore if the adjoint operator \P := $T = 
J _1 M of this hereditary operator $ = MJ -1 maps a gradient vector /o = to another 
gradient vector f\ = tyfo = then all vectors ^"/o, n > 0, are gradient, i.e. there exist 

functionals H u i > 2, such that & f Q = i > 2 (see [§ @ [§ for more information). 
Then it follows from bi-Hamiltonian structure that all systems of evolution equations in 
the hierarchy m = & n fo, n > 0, commute with each other, i.e. [$ m /o, 3> n /o] =0, m, n > 0, 
and they have infinitely many common conserved densities being commutative under two 
Poisson brackets. Therefore Hamiltonian pairs can pave a way for constructing integrable 
systems. However the problem still exits and just turn to how to find Hamiltonian pairs. 

In this paper, what we want to develop is to propose a possible way to generate Hamil- 
tonian pairs. We are successful in constructing bi-Hamiltonian coupled KdV systems in 
such a way. The paper is organized as follows. We first concentrate on the techniques for 
extending Hamiltonian operators from lower order to higher order of matrix, motivated 
by the idea in @ |]. Then we will go on to analyze a class of Hamiltonian pairs which 
can yield hereditary operators and eventually new bi-Hamiltonian coupled KdV systems. 
A few of concluding remarks are given in the final section. 
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2 Extending Hamiltonian operators 



Let us specify our dependent variables 

u k = {ul,---,ul) T , l<k<N, 

I T T \T i 1 q l q \T 

and introduce a condition 

J' k {u k ) = J[{u l ), l<k,l<N, (2.1) 

on a set of given Hamiltonian operators Jk{ u k) '• A q (u k ) —> A q (uk), 1 < k < N. This 
condition requires a kind of linearity property of the involved operators with regard to their 
dependent variables. Such sets of Hamiltonian operators Jk{uk) do exist. For example, we 
can choose 

Jk(u k ) = ~ld 3 + 2du k d- 1 + 2u k , d = ^-, 1 < k < N. 
4 ox 

The problem that we want to handle here is how to generate Hamiltonian operators 
starting from a given set of Jk(u k ), 1 < k < N. A simplest solution is to make a big 
operator to be J(u) = diag(Ji (u% ),■■-, Jn(un))- What we want to develop below is to 
propose more general structure of Hamiltonian operators. To the end, we introduce the 
following structure of candidates for Hamiltonian operators 

N 

^) = (E4 J ^)) W (2-2) 
fc=l 

where {c^ | i,j, k = 1, 2, • • • N} is a set of given constants. Obviously this big operator J(u) 
may be viewed clS cl linear operator 

J(u) : A Nq (u) = A q {u) x • • • x A q (u) A Nq {u) = A q {u) x • • • x A q {u) , 

v ' V v ' 

N N 

where a vector function of A q {u) depends on all the dependent variables u\, ■ ■ ■ ,un, not 
just certain dependent variable u k , by defining 

rp N 

Ja = (((Ja)i) T , • • • , ((Ja) N f) , (Ja) t = ^ c^J k (u k )a j} l<i<N, (2.3) 

j,k=i 

where a = (af , ■ ■ ■ , aJj) T , ol% G A q (u), 1 < i < N. To guarantee the skew-symmetric 



property of the big operator J, a symmetric condition on {cfA 



c% = c%, l<i,j,k<N (2.4) 

suffices. The following theorem provides us with a sufficient condition for keeping the 
Jacobi identity (|Q|), 

Theorem 2.1 // a/Z J k (uk) '■ A q {u k ) — > A q (u k ), 1 < k < N, are Hamiltonian operators 
having the linearity condition ( \2.1[ ) and the constants c^- ; 1 < i,j,k < N, satisfy the 
symmetric condition ( OO j and the following coupled condition 

N N 

X) 4 C W = E 4^-, 1 < i, 3, hn<N, (2.5) 

k=l k=l 
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then the operator J(u) : A Nq {u) -> A Nq (u) defined by (O) and ( O ) is a Hamiltonian 
operator. 

Proof: We only need to prove that J(u) satisfies the Jacobi identity (|Q|), because the 
linearity property and the skew-symmetric property of J have already been shown. Noting 
that A q {u) is composed of column function vectors, we suppose for a, 0,7 G A Nq {u) that 

a = (al,---,a%) T , (3 = (/f , • • • , /?£) T , 7 = ( 7l T , • • • , 7 £) T , ai,A,7i G A q (u), 1 <i < N. 

Moreover we will utilize a convenient notation (X)i = Xi, X{ G A 9 (u), 1 < i < iV, when 
a function vector X G ^4^ 9 (u) itself is complicated. 

First from the definition fl2.2|) , we have 

iV 

NxN 



1 V 

J'(n)[J/5] = (^44K)[(J/3) A 



.*J 

fc=i 

A 7 " 

(J'(n)[J/5] 7 ) l = £ 44(n fe )[(J/5),] 7 „ l<i<N. 
j,k=i 



Then taking advantage of the concrete definition ( |2.3| ), the linearity condition ([O]) and 
the coupled condition fl2,5|) , we can make the following computation 



< a, J'{u)\jp\y > +cycle(a, /?, 7) 

JV 

£ c*- < «i, 4(ii fe )[(J/3) fe ]7j > +cycle(a,/3, 7 ) 
ij',fe=l 

N N 

c ij < a>i,J' k (u k )[ Cki J n( u n)Pi]lj > +cycle(a, /3, 7) 

i,,7,fe=l (,n=l 

12 c ij c fcj < a *> 4( u fe)[ J n(«n)A]7j > +cycle(a, (3, 7) 

i,j,k,l,n=l 
N N 

12 £( c ?j c w) < a *> 4( n fe)[ J n(^n.)A]7j > +cycle(a,/3,7) 

i,j,l,n=l k=l 
N N 

12 Y.(°ij C kl) < a h Jn(Un)[Jn(u n )Plhj > +CVcle(a, /?, 7) 
i,j,l,n=l k=l 
N N 

12 12( c ij c kl) < «i> ■^(«n)[J r n(«n)A]7j > 
i,j,l,n=l k=l 

N N 

+ J2 12( C ij C kl) < Pi,J' n {Un){Jn{u n )ll]0Lj > 
i,j,l,n=l k=l 
N N 

+ 12 Y,( c iJ c ki) < 7i, ^K)[JnK)a/]/3j > 

i,j,l,n=l k=l 
N N 

12 12( c ij c ki)(< a h J' n ( u n)[Jn(u n )Pihj > +cycle(a i; (3 t , 7j 

i,j,l,n=l k=l 
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Recall each Jiiui) is Hamiltonian and it follows from the above equality that the big 



operator J defined by ( |2.2[ ) and ( |2.S| ) satisfies the Jacobi identity (1.3). The proof is 
completed. I 

We remark that if we consider the coefficients {c^ } to be the structural constants of a 
finite-dimensional algebra with a basis ei, e2, • • • ,ejv as follows 



N 



(2.6) 



fc=i 

then the symmetric condition ( |2.4[ ) and the coupled condition ( |2.5[ ) can become 

ej * ej = ej * 1 < i, j < iV, (2.7) 

(ei * e,) *e k = (e k * e { ) * e j: l<i,j,k< N, (2.8) 
respectively. They reflect two specific properties of the related algebra. 

Apparently basic scalar Hamiltonian operators satisfying the linearity condition (|2.l| ) 
can be the following set 

Ji( Ui ) = ad 3 + did + 2u ix + 4uid, l<i<N, (2.9) 



where d = d/dx and q, dj, 1 < i < N, are arbitrary constants. Of course, matrix Hamil- 
tonian operators having the linearity condition Q2.1| ) may be chosen. Actually, such sets of 
Hamiltonian operators may be presented directly from the above operators by Theorem 
2.1 or by perturbation around solutions as in Refs. M |T( 



In what follows, we example applications of Theorem 2.1 to two specific choices of 



{ck}. The analysis below is quite similar to the one made for the extension of hereditary 



operators in Ref. 

Example 1: Let us first choose 



Cjj Oi+j,k—p: 



(2.10) 



where p is an integer and 5 k i denotes the Kronecker symbol. The corresponding big operator 

J p+3 (u p+ 3) 



formed by (2.2) becomes 

Jp+2(Up+2) 



J(u) 



J p+3 (n p+ 3) 



Jp 

Jp+N+2{u p +N+2) 
Jp+2N(Up+2N) 



(2.11) 



Jp+N+l(u p +N+l) Jp+N+2{u p +N+2) 

where we accept that Ji(ui) = if i < or i > N + 1. 

For two cases of —2N +1 < p < —N and — 1 < p < N — 2, the coupled condition (|2,5| ) 
can be satisfied, because we have 



TV 



N 



k=i 



k=i 



i, 

0, 



when n — i 
otherwise. 



3-1 = 2p, 
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While proving the above equality, we should keep in mind that we have 

1 <i+j+p = n — I— p< N, 1 <i+l+p=n— j — p < N, 

when n — i — j — I = 2p. But for the case of —N < p < —1, upon choosing i = —p + 1, j = 
—p — 1, n = I = 1, we obtain 

N N 

^ c ij c kl = h c li c kj = 0' 
k=l k=l 

and thus the coupled condition (|2.5|) can not be satisfied. 

Note that when p > N — lorp< —2N, the resulting operators are all zero operators. 
Therefore among the operators defined by ( |2.11 ), we can obtain only two sets of candidates 
for Hamiltonian operators 



J(u) 







Jp+N+l(u p+ N + i) 
Jp+N+l( u p+N+l) Jp+N+2{ u p+N+2) 



. Jp+N+l{Up+N+l) Jp+N+2{Up+N+2) 



J(u) 



Jp+2N{Up+2N) 

-2N + l<p<-N, (2.12) 

Jp+2{u p +2) Jp+?,{Up+z) ■•• Jp+N+l(Up+N+l) 



Jp+3(u p +3) 



Jp+N+l(u p +N+l) 



-1 <p < N - 2, 



(2.13) 



where we still accept that <&i(ui) = if i < or i > N+l. These two sets of operators can be 
changed to each other by a simple transformation (ui, U2, ■ ■ ■ , ujy) «-> (un, wjv-Ij " " * > '"l)- 



Example 2: Let us second choose 



hi, l = i + j-p (modiV), 



(2.14) 



where 1 < p < N is fixed and 5ki denotes the Kronecker symbol again. In this case, we 
have 



^ n _A kn _ f !> wheni + j+/-n = 2p (mod AT), 

/ y CijCkl / j ^li^kj 



(2.15) 



fe=i fc=i I 0, otherwise, 

which allows us to conclude that the coupled condition (gTa) holds, indeed. Therefore we 
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obtain a set of candidates for Hamiltonian operators 

J2-p(u2-p) J3-p{u3-p) 
J3-p(u3-p) 



J{u) 



Jn-p+i{uN-p+1) JN-p+2{uN-p+2) 



Jn-p+i(un^ p+ i) 

Jn- P +2{uN-p+2) 
J2N-p(u2N~p) 



(2.16) 



where we accept Jj(uj) = Jj(v,j) if i = j (modiV), while determining the operators in- 
volved, for example, J 2 - P (u 2 - P ) = Jn(un) when p = 2. A special choice of p = 1 leads to 
a candidate of Hamiltonian operators 



J(«) 



J2(u 2 ) 



Jn(un) 
Ji(ui) 



(2.17) 



Jn(un) Ji{ui) ■■■ Jn-i(un-i) 

This operator will be our starting point for constructing new integrable systems in next 
section. 



3 New bi-Hamiltonian coupled KdV systems 



From now on, we focus on the candidate of Hamiltonian operators defined by (|2.17| ). Let 
us pick out the Hamiltonian operators formed by (|2.17|) under a choice of (|2.9|). Then the 
following Hamiltonian pair can be engendered 



J(u) = Ad 



(12 



a N 
ai 



0. 



(3.1) 



M(u) 



aN a\ ■ ■ ■ ajv-1 
Mi(«i) M 2 (u 2 ) ••• M N (u N ) 



M 2 (u 2 ) 



Mi (ui) 



M N (u N ) Afi(«i) ••• M7v-i(njv_i) 



(3.2) 



where d = d/dx, u = (ui, u 2 , ■ ■ ■ , «at) t , aj = const., 1 < i < N, and the operators 
Mj('Uj), 1 < i < iV, are given by 



Mj(uj) = Cj(9 3 + did + 2^ + 4uj9, Ci,d{ = const., 1 < i < N, 



(3.3) 



which are all Hamiltonian. 
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We assume that the constant matrix A is invertible to guarantee the invertibility of J, 
and its inverse matrix is given by 



B 



bi b 2 
h 



b N 



(3.4) 



b N bi ■■■ b N -! 
where bi, 1 < i < N, can be determined by solving a specific linear system 

a-ih + a 2 b 2 H h a N bN = 1, 

(i2h + a 3 b 2 H h a N bN-i + ai6v = 0, 



ffljv^i + a\b 2 H h a N _ib N = 0. 

Now the resulting hereditary operator = MJ~ l reads as 



Mi(«i) M 2 (u 2 ) 
M 2 (u 2 ) 

M2v(«jv) Mi(ui) 

$i(iti) $ 2 (u 2 ) 

$2(^2) 



M N (u N ) 
Mi(«i) 

M/V-lOujV-l) 
$>Ar(UAr) 

$l(«l) 



&1 
&2 



61 6 2 
b 2 

b N b 1 
b 2 



b N 
h 

bN-i 
b N 



d- 1 



bN h ■■■ bN- 



(3.5) 



with 



§i{ui) = Miiu^d' 1 = ad 2 + di + 2u ix d~ 1 + Au h 1 <i < N. (3.6) 
This hereditary operator can be rewritten as a concise form 



<J?(u) = MJ- 1 = (^ph-i+^kiuk)^ , 



(3.7) 



where 6j = bj if i = j (mod N) . It is also an example satisfying the extension scheme of 
hereditary operators in Ref. [pi]], because upon setting c^- = bk-i+j we have three equal 
sums for all 1 < i, j, l,n < N: 



N 



N 



N 



/ 1 ^ij^kn 


= ^ bk-i+jbl-k+n 


— ^ t bmbl+i—j+n—mi 


k=l 


k=l 


m=l 


N 


N 


N 


C ik c kj 


= ^ bi-i + kb n -k+j 


= ^ t bmbl+i—j+n—mi 


k=l 


k=l 


m=l 


N 


N 


N 


/ , ^in'-kj 
k=l 


= ^ bk-i+nbl-k+j 
k=l 


— ^ ' bmbl+i—j+n—mi 
m=l 
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which are sufficient for <3?(-u) to be hereditary (see 

We now turn our attention to investigating the nonlinear systems which can be resulted 
from the above Hamiltonian pairs. Such first system can be the following 



ut = &u x = MJ 1 u x , 

which can be represented as 

TV 

uu = ^2 b k ^ i+j (c k u jxxx + d k u jx + 2u kx uj + Au k u jx ), 1 < i < N, 
k,j=i 

where we again accept hi = bj if i = j (mod N) . It is easy to find that 



(3.8) 



(3.9) 



fo := J u x = Bu 



5H n 



1 



^ , H = j H Q dx, H = -u Bu. 



(3.10) 



Go ahead to check whether or not the next vector defined by 

^ 2 {u 2 ) ••• *n(u n ) 
V 2 {u 2 ) '•• 



h ■■= * /o = * + /o = B 



^ N (u N ) *l(ui) ••• tfjV-l(«iV-l) 



Bu 



(3.11) 



is a gradient vector, where ^(ttj) = $|("Uj) = d 1 Mj(ui) , 1 < i < N. That is true, indeed. 
Actually we have 



fi = */o = -j^, Hi = J H 1( ix, 

e 1 ( Ul ) e 2 (u 2 ) 
e 2 (n 2 ) 



H x = (Bu)" J 



9i(ui) 



Bu 



N N ^ 

52 (52 b i+l-ih+j-i)ui(-c k Uj xx + ~d k Uj + U k Uj), 

j,k,l=l i=l 



(3.12) 



where the operators 0j(uj), 1 < i < N, are given by 



1 12 2 

®i{ui) = -Cid 2 + -di + -Ui + -d~ 1 u i d, l<i<N, 



(3.13) 



and hi = bj if i = j (mod iV). We can also choose the energy form for the functional H\: 



N N 



f 1 1 

H i = J Hidx, Hi = 22 (52 bi + l ~ lbk +j~^~2 CkUlxUjx + 2 dkUlUj + UlUkU i)- ( 3 - 14 ) 



j,k,l=l i=l 
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Therefore according to the Magri scheme H 01 PI; there exist other functionals 
H n , n > 2, such that ^f n fo = t^S n > 2. All such functionals can be generated by 
computing the following integrals 

H n = f f 1 < (y n f )(\u),u >d\dx, n > 0. (3.15) 



o 



Further we can obtain a hierarchy of bi-Hamiltonian equations 

Ut = Kn := (Hu)Y +l u x = J^±i = M^, n > 0, (3.16) 

ou ou 



which includes the nonlinear system fl3.9|) as the first member. It follows that they have 
infinitely many commutative symmetries {-ftT m }o° and infinitely many commutative con- 
served densities {i^ m }o°. All systems of evolution equations can reduce to KdV equations 
once Uj = Ci = dj = 0, j ^ 1, are selected. Therefore they are all iV-component coupled 
KdV systems. 

Let us now work out a concrete example for a choice of A with a± = 1, Oj = 0, 2 < « < 

N. 

In this case, we have b\ = 1, = 0, 2 < i < N. Then the first Hamiltonian structure 
is given by 

dtut = 8 X ^, d t u k = cL 6Hl , 8 X = 2 < k < N. (3.17) 

OUi 0U N+ 2-k OX 

This first Hamiltonian structure has the momentum: 

1 f N 

P = H = - [uj + Y^ u k u N+2 -k]dx, (3.18) 

Z J k=2 

which is of quadratic form with respect to its N Casimirs (annigilators of the first Poisson 
bracket) F k = J u k dx, k = 1, • • • , N. The conservation law of momentum is 

N rfV 

d t H 1 = d x [J2u klr L -F}, (3.19) 

where H\ = j H\dx, d x F = J2k=i J^ u kx- For simplicity in the selection ( |3.6| ) we can put 
di = 0, because we can eliminate those constants in the expressions of <3?j by making shifts 
Ui — > Ui — di/A. At this moment our coupled KdV system reads as 

m N m N 

dtu m = dx[dx(^2 c m+i-kUk+ Cjy+m+i-fc^fc)+3 u k u m+ i_ k +3 u k UN+ m +i- k ], 

fc=l fc=m+l fc=l k=m+l 

(3.20) 

which has the following Hamiltonian for the first Hamiltonian structure 

Hi = [ H x dx (3.21) 



^ N 2V-1 N+l-m 

Hi = --[ciwj+WiJ2c k W N+ 2-k+ C m ^ 

k=2 m=l k=2 
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N N N 

+ ^2 Cm W k W 2 N-m+3~k] + uf + 3ui ^ U k U N+2 -k 

rn=3 k=N+3-m k=2 
N-l N+l-k N-3 m 

+ Y Uk X! u mUN+3-k-m + Y U m+3 Y 

^N—k^N+k—mi 

k=2 m=2 m=0 k=0 

where w k = d x u k , 1 < k < N, and N is assumed to be greater than two for producing 
nontrivial systems. The second Hamiltonian structure can be determined by its recursion 
operator 

' 5H N 5H 

d tUl = T $ i+1 - k d x —^ + T $ N+1+i - k d x —?-, 1 < i < N, (3.22) 



where <3?j = Cid% + 2(2uj + Wid x 1 ), 1 < i < N. The second Hamiltonian structure has the 
momentum F[ = J u\dx and the Hamiltonian Hq. 

The relationship between the gradients of H k and Hk+l determined by recursion oper- 
ator yields a possibility for constructing an infinite set of conservation laws and commuting 
flows by iterations. At each step we need to compute integrals to construct Hk+l- They 
can be done by using the formula ( |3.15 ) in variational analysis. 

Moreover we have an alternative way to construct an infinite set of conservation laws 
and commuting flows. Let us introduce an eigenfunction problem for the recursion operator 

[*(«) - X]u T = 0, (3.23) 

i N 
i.e. Y ^i+i-kdrUk + ®N+i+i-kd T u k = Xd T u i: 1 < i < N. (3.24) 

fe=l k=i+l 

This eigenfunction problem can be equivalently rewritten as 
f d T u k = d x v k , 1 < k < N, 

i N 

Y ®i+i-kd x Vk + Y ®N+i+i-kd x v k = Xd x Vi, 1 < i < N. 

k=l k=i+l 

If we use formal series near A — > oo 

d T = d t0 + jd tl + -^d t2 + ■ ■ ■ , (3.25) 
Vi = ^° ) + ^ 1) + ^f ) + -, (3.26) 

we can obtain directly from the above eigenfunction problem an infinite set of commuting 
flows. Furthermore we can use formal series near A — > 0. Then the first resulting nontrivial 
commuting flow is a coupled long wave equation 

i N 

$i+i-kd T u k + $N+i + i-kd T u k = 0, 1 < i < N. (3.27) 

k=l k=i+l 



This is an N-component generalization of long wave equation (see [12] for more information 
about long wave equation), which commutes with the KdV equation. 
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4 Concluding remarks 



We have proposed new coupled KdV systems possessing bi-Hamiltonian structures by 
extending Hamiltonian operators from lower order to higher order of matrix. Clearly we 



may make other choices of Ji(iii) to give more results based on Theorem 2.1 
Compared to the well-known coupled KdV systems (for example, see [O]) 



Ut = <& n u x , $(«) 



■■• $i 

1 ••• $ 2 



1 $ 



N 



n > 0, 



(4.1) 



and quite new couple KdV systems introduced in (14] 



u t = <S> n u x , $(«) 



b N ®l 

&jv-i$i + b N <$> 2 b N <5>i 



bi$i H h b N <Z> 



A' 



n > 0, 



where u = (u\, U2, • • • , un) , &i, 1 < i < iV, 
$ i = ^(ui), 1 < i < N, are still defined by ( |3^ ). 
similar nice integrable properties, for example, 



6iv_i$i + b N <S>2 b N $i 

(4.2) 

are arbitrary constants except bjy / 0, 
our new coupled KdV systems ( 3.16| ) have 
bi-Hamiltonian structures, dispersionless 



limits having bi-Hamiltonian structures (see [15| for the case of the well-known coupled 
KdV systems). But there exist differences among the structures of recursion operators 
corresponding to these three hierarchies of coupled KdV systems. The bi-Hamiltonian 



structure (see [|T4|) of the hierarchy (4.2) can similarly be derived from the Hamiltonian 
operators formed by (|2,12j ) in the first example of the second section. This is why we didn't 
deliver above a detailed analysis for constructing integrable systems starting from ( 2.12j ) 
or equivalently from ( |2.13| ). However we don't know whether or not two new coupled KdV 
hierarchies (|4.2| ) and ( |3.16| ) have other integrable properties, for example, Lax pairs like 

yD- 



In terms of the existence of recursion operators, other integrable couple KdV systems, 
say, Jordan KdV systems, have also been derived (see for example [16] [0 ][il| |l9|]). A nat- 
ural question is whether there exist other hierarchies possessing bi-Hamiltonian structures 
among those coupled KdV systems. This will enrich the content of Hamiltonian theory for 
coupled KdV systems. 
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